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SUMMARY 

The Orbiting Geophysical Observatory (OGO) is an actively 
oriented satellite -which will carry up to fifty experiments. The 
experiments will be oriented relative to the Sun, Earth or the orbit 
plane. This document presents methods of analysis for computing orbital 
data, spacecraft angle data and spacecraft heat input data for OGO, and 
gives specific results for the S-l+9> Eccentric Orbiting Geophysi-cal 
Observatory (EGO). The orbital data include the transient behavior 
of true anomaly, the flight path angle, the distance from the center 
of the Earth to the satellite, the angle between the line of apsides 
and the ecliptic plane, and the angle between the projection of the 
line of apsides on the ecliptic plane and the Earth-Sun line. The 
spacecraft angle data include the transient behavior of solar array 
angle, the orbital plane experiment package (OPEP) angle and the OPEP 
velocity vector angle. The transient heat input data include Earth 
emitted heat, direct solar heat and Earth reflected solar heat inputs 
to all the faces of the mein OGO box, the solar array and the OPEP. 


EGO ORBITAL PARAMETERS 
AND HEAT INPUTS 

by 


B. E. Montgomery 
S. J. Paddack 
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Goddard Space Flight Center 
Greenbelt, Maryland 

I. I NTRODUCTION 

This report presents orbital data, spacecraft angle data and 
heat input data for EGO. A detailed description of the methods used 
in the calculations is also given. The calculated results are con- 
sistent with the "EGO Launch Window Study", NASA Report Number X6h}-62-225 
(Reference l), insofar as the perigee restraint is concerned. A slightly 
different set of initial conditions was used in Reference 1, but the 
launch window map will still give a good estimate of launch times 
consistent with the perigee restraint. The data presented in this report 
are for launch times for which the perigee will not go below its initial 
value for a whole year. The injection parameters for this report were 
taken from Reference 2. 



Azimuth 

north 


the compass heading measured clockwise from 


The semi-major axis of the satellite’s orbit. 

1 A.U. = 1.496 x 10 B kilometers = the semi -major axis 
of the Sun’s orbit relative to the Earth (Reference 6). 

(i^ x i e ) • by definition. 

The number of days from 0.0 hours Universal Time on 
January 1, 1961 to the time in question. 

An elemental area of the Earth's surface. (See Figure 3 ). 

The satellite’s eceentric anomaly. 

The eccentricity of the satellite’s orbit. 


0 . 0l£>7259 » the eccentricity of the Sun's orbit relative 

to the Earth. (Page 50 of Reference 7 ) . 


oS 


— - = the reflected solar radiation factor. 


A 


— = the Earth emitted radiation factor. 


% 

S 


n • r, 


sv 


l r svl 


the direct solar factor. 


= The integrand of the integral for q a . (See page 26 ) 


NOMENCLATURE (Coot) 


p) = The integrand of the integral for q^,. (See page £5 .) 

G Mg = 3.986032 (±0. 000030) x 10 s km 3 , see” 2 = 

the gravitational constant of the Earth. ‘ 

(Reference 6). 

g* « - Mean anomaly of the Sun relative to the Earth. 

_ X VE 

= The angular momentum vector. 

the distance from the reference sphere 
(of radius equal to the Earth's equatorial radius) 
to the spacecraft measured along the geocentric 
radius vector. 

i = Inclination = the angle between the satellite's orbit 

, plane and the equatorial plane. 

i, ,5, k * The unit vectors which correspond to X, Y and Z, respectively. 

i^, k^ - = The unit vectors which correspond to x^, y fe , and z fe , 

respectively. (See Figure l). 

* The unit vectors which correspond to the x . y , and z 

. e' e e 

coordinates, respectively. (See Figure l). 

ip, «}p, kp =« The unit vectors which correspond to y^, and z , 

respectively., (See Figure l). 



n 


H 


m 


r VE X r VE 


Height = 
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NOMENCLATURE (Cont) 


k* 


L* 


p, Q, R 

P Y» p z 

St* Sc* ^z 


The unit vector vhich is normal to the ecliptic plane 
and points north. 

u)» + g* = the mean longitude of the Sun. 

The unit outward no rmal of the surface in question. 

The unit vectors vhich correspond to x^, y^, and z^, 
respectively. (See Figure 2). 

The components of the P unit vector in the (X, Y, Z) 
system. 

The components of the Q unit vector in the (X, Y, Z) 
system. 



= The heat flux, incident upon the flat surface, which originates 
at the Sun and reflects from the Earth to the vehicle. 




The heat flux, incident upon the flat surface, which radiates 
from the Earth. 


% 


The heat flux, incident upon the flat surface in question, 
whose source is the Sun. 


The equatorial radius of the Earth. 

The components of the R unit vector in the (X, Y, Z) 
system. 


V ^ 


hi* 


IiOI-ERCLATUKE (Cont) 


SE 


Tiie position vector which points from the center of the 
Earth to the Sun. 



VE 


l VS 


The position vector from the Earth to the vehicle. 

The velocity of the satellite relative to the Earth. 

The position vector which points from the Sun to -the 
vehicle. 


S 


U SE 


= I.97 ± 0.01 cal cm" 2 minutes 1 = 0.137^ watt, cm 2 

= Flux of total radiation received outside Earth's atmosphere 
per unit area at mean Sun-Earth distance. (Reference 12) 

The unit vector which points from the center 
of the Earth to the Sun. 



v = The speed of the vehicle. 

S ' 


X, Y, Z = Cartesian coordinates with the origin at the center of 

the Earth, the X-axis in the direction of the vernal 
equinox, the Z-axis along the north pole of the Earth, 
and the Y-axis forms a right hand set. 


This coordinate system is fixed to the main 0G0 box. 
(Figure l). 

x . y . z = The OFEP coordinates and are defined with respect to 

c e 

the (x fe , y fa , z b ) system by the OPEP angle * e . 
(Figure l). 


V y b- 2 b 
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The paddle coordinates and are defined with respect 
to the (x^, z^) system by the paddle angle cp . 

(Figure l). 

Cartesian coordinates with the origin at the center of 
the Earth, the x^. - axis is in the direction of perigee, 
the y^ - axis is advanced 90° in the direction of motion, 
and the 2 - axis is in the direction of the satellite’s 
angular momentum vector. (Figure 2). 

The albedo » the fraction of the solar constant S which 
is reflected away from the Earth. The value for a is 
about 0.54. (Reference 12). 


The angle between the vector u<,g which points from the 
center of the Earth to the Sun and the outward-directed 
normal of the Earth’s elemental surface area cLA. 

(Figure 3). - 


Flight path angle = the angle between the velocity vector 
and the plane normal to the radius vector passing through 
the vehicle. 




'- 9 O 0 by definition. 


The angle between the equatorial plane and the ecliptic 
piane. * 

The angle between the outward-directed normal of dA and 
the vector 7 • (Figure 3). 



NC0SHCLATUBE (Cont) 


6 


6 e 


% 

A 


X 


X’ 




v 

V* 


= . Terrestrial longitude * the angular distance from the 

Greenwich meridian, measured eastward along the equator 
to the meridian plane of the vehicle. 


= OPEP - velocity vector angle * the angle between the 

vectors T and T xre , . 

B VJb 



, (Figure 3). 


= 0.02215 watts cm 2 = the emittence per unit area of 

the surface of the Earth. 




cos 


(-9 


VE 


KeI 


• n'y - the angle between the vector 


- r^j, and the unit vector n . (Figure 3) 


as v* + U) 1 • True longitude of the Sun = the angle between 

the positive X-axis (Vernal Equinox) and the vector which 
points from the Earth to the Sun . It is measured 
in the direction of motion from the X-axis. 


= The angle between the outward-directed normal of the plate 
n and the vector (-^) . (Figure 3). 

« The satellite's true anomaly. 

= True anomaly of the Sun * the angle between the vector 
which points from the center of the Earth to the Sun's 
perigee and the vector which points from the center 
of the Earth to the Sun. It is measured in the direction 
of motion of the Sun. 



The angle between the projection of the (-P) on the 
ecliptic plane and . 

The vector which points from the Earth 's elemental area 
dA to the vehicle. (See Figure 3). 

i 2 

5.735 x 10“ watt cm”" 2 °K ~ 4 = the Stephan- Eoltzman 

constant. 

Time from the injection into orbit. 

Geocentric latitude = the angle at the center of the 
Earth between the radius through a given point and the 
equatorial plane. 

The angle between the plane formed by the vector r VE 
and Tgg ; and the plane formed by the vector r^ 
and n. (Figure 8). 

sin f 1 ) = the maximum look angle for which 

v IT y 

the satellite can see the surface of the Earth. 

(Figure 4). 

The solar array angle * the angle measured counter 
clockwise from the vector to the j vector when 
observed from the positive end of the vector. 
(Figure l). 

The minimum angle between the vector pointing from 
the center of the Earth to apogee and the ecliptic 
plane. 


NOMENCLATURE (Cont) 


Q 


0) 


«* 


r 


**> 


= The OEEP angle = the angle measured counter clockwise 

from the vector i, to the vector i , when observed 
. d e _ 

from the positive end of the vector . 

= Longitude of the ascending node = the angular distance 

from the vernal equinox measured eastward in the equator- 
ial plane to the point of intersection of the orbit 
plane where the satellite crosses from south to north. 

= The argument of perigee = the angular distance measured in 
the orbital plane from the line of nodes to the line of 
apsides. 

= The argument of perigee of the Sun « the angle between 
the positive X-axis (Vernal Equinox) and the vector 
from the center of the Earth to the Sun*s perigee 
measured in the direction of motion from the X-axis. 

= The tiue rate of change of f> 

= The tiue rate of change of uj • 
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III. ASSUMPTIONS 


A. INITIAL CONDITIONS 

The initial conditions are assumed to be the same as the 
Agena burnout conditions (Reference 2 ) namely: 


ft 

- 

-20.280 

degrees 


geocentric latitude 

e 

as 

111.902 

degrees 

= 

terrestrial longitude 

h 

s 

281621*6 

meters 

. m 

geocentric height 

V 

s 

= 

1064p.27 

meters per second ** the speed 
of the vehicle 

A c 

* 

66055 

degrees 

- 

azimuth 

y 

S 55 

1 . 1+7148 

degrees 

* 

flight path angle 


B. RESTRAINTS ON LAUNCH TIME 

For the purpose of generating results for this report the 
only restr ain t on launch time is the one on perigee latitude. This 
restraint states that the perigee altitude will not go below its initial 
value for one year. The allowable launch times were chosen from the 
launch window map in Reference 1 . A slightly different set of initial 
conditions was used in Reference 1, but the launch window map gives a 
good estimate of launch times consistent with the perigee restraint. 

IV. MATHEMATICAL MODEL 

A. IDEAL CONTROL LAWS AND DEFINITION OF COORDINATE SYSTEMS 
Figure 1 presents the body coordinate systems. The 
V coordinate system is fixed to the main box. The solar 

paddle onordi nnte svstem x . v . z is defined with respect to the 

* ~ ~ p' -p' p 
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x b r y b , 2 b system by the paddle angle t' p . The y axis is the outward- 
directed normal to the surface of the paddle which bears the solar 
cells. The OPEP (Orbit Plane Experiment Package) coordinate system 
x e* y e' z e is defined with respect to the x^ y , z fe system by the 
OPEP angle The x axis is the outward-directed normal of the 

surface of the OPEP which is facing in the direction of motion. 

It is assumed that the CGC spacecrafts obey ideal control laws. 

A definition of these control lavs fo llo ws: 

1. The 4 - face will be aligned with radius vector from 

the center of the Earth to the satellite. (Page 23 of Reference 3). 

In equation form this states that 



2. The solar array (x^-z^ plane) will be normal to the vector 
from the Sun to the satellite. In equation form this states 
that 
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3. The Sun will not shlr.e or the + y - face of the raair. box 
(Page 20 of Reference 4). In equation form this states that 

r VS * "*b = r VS ^ \ X /S ^ = 0 

i+. The positive x_ - face of the orbit p lane experiment 
package (OFEP) shall look forward in the plane of the orbit. 
(Page 24 of reference 3). In equation form this states that 


and 



i 

e 


* 0 


i 

e 


H m 

P „ * 5,1 


Where H = r._ x r,^ is the angular momentum vector. With these 

four control lavs three orthogonal coordinate systems may be 

established; box (I b , i^), solar array ( l v ^,/'yand OPEP 

(I . 3 , k ) . The vectors k, and 3 are outainea directly from 
e* e e d p 

laws 1 and 2 respectively. Since T b is perpendicular to both 

T, and it follows that • 
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The sign of the expression for i. is determined from control 

D 

law number The vector j is determined from 7 - K x T, . 

o bob 

By construction k = k . Also, by construction i - i . Then it 
ah p " b 

follows that k = j x i - The vector i is determined by control 
P p P _ e 

law number h r and lastly, - k^ x i e * 

B. DEVELOPMENT OF POSITION AND VELOCITY VECTORS. 

1. Earth-Sun Vector r 

oil* — 

This is the vector which points from the center of the 
Earth to the center of the Sun. 

r SE = r SE U SE 

• where r^-, is the magnitude of r , and u^g is the unit vector 
which points from the center of the Earth to the center of 
the Sion. The value of r ^ may be calculated with the 
expression (page l6h of Reference 5) 

a* (l-e* 2 ) 
r SE 1 + e’ cos n* 



where 


a’ = 1 a.u. = 1.^9599 x 10 8 kilometers = the semi-major, 

axis of the Sun’s orbit relative to the Earth. (Reference 6). 

e’ = O.OI67259 * the eccentricity of the Sur's oibit relative 

to the Earth, (page 50 of Reference 7). 

v* = the true anomaly of the Sun relative to the Earth 
^g* + 2 e’ sin g' , since e’ is small, (page 171 of 
Reference 5) . 

g* = - 3.2419100 x 10~ 2 + 1.7201970 x 10" 2 d = the mean 

anomaly of the Sun relative to the Earth. (Page 498 of Reference 7). 

d = the number of days from 0.0 hours Universal Time on 

1 January, I96I. (Page 498 of Reference 7). 

* 

The unit vector Ugg is given by 

Ugg « i cos X' + j cos c sin X* + k sin e sin X* 
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where 


e = The angle between the equatorial plane and the 
ecliptic plane. 

X’ = True longitude of the Sun. 

= uj' + '*’ = u)‘ + g* + 2e sin g* = L* + 2e sin g* 

« k.Eeh 117 + 0.017202791 a 

+ 0.033^5099 sin P -0.032^191 4 0.01720197 d J 
<#' = Argument of perigee 

L* = 4 g* = Mean longitude of the Sun. 

A more convenient form for the -unit vector Ug^ is 


~ p * Q) Q - (Ugj, * R) 


R 


"SB 


= (P v cos ) T + P v cos e sin X* + p sin e sin X*) P 


+ cos ).' 4 Qj cos e sin x* 4 ^ sin e sin X') Q 


+ (K, cos X« 4 R^. cos e sin X* 4 sin e sin X*) R 
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where 


of the unit vector P 
of the unit vector Q 
of the unit vector R 

2. Earth -Vehicle Vector r^ 

The vector r.^, points from the center of the Earth to the vehicle 
The Launch Window Program cannot be used to calculate the vector r^ 
as a function of time for tracking purposes. It can, however, give 
an adequate estimate of the vector r^ relative to other vectors 

r SE* r SV J i b J k b' 1 p’ k p* 1 e > k e* etC 

as a function of time from perigee. 

It is assumed that at any particular time t for which the Launch 
Window Program has calculated the orbital elements (a, e, i, ob, 0), 
the satellite is at perigee. At times subsequent to t, up to one 
orbital period, two body motion is assumed. The position vector may 
be written in the form (Pages 5^ and 65 of Reference 8) 

Ta (cos E - e) 4 Qe ^ 1-e 2 sin E 



P^., P^, are the components 


St' S' S; are the components 


R x» r y j are the components 
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where, 

a = The serai -major axis cf the satellite’s orbit 

E = The eccentric anomaly ox" the satellite’s orbit 

e = The eccentricity of the sate lli te’s orbit 

P = The unit vector which points from the center of the 
Earth to the satellite’s perigee 

Q = The unit vector which is rotated 90° in the plane of the 
satellite’s orbit in the direction of motion. 

These assumptions can be made since the relative locations of 
the Sun and the Earth with respect to the orbit do not change enough 
in one EGO period to be seen as a first order effect on any of the 
parameters presented in this report. 

The errors in orbit shape and orientation are allowed to accumulate 
for only one orbital period, since the orbit is rectified after each 
period. The Sun moves approximately one degree per day relative to 
the Earth; the orbit plane precesses about 0.65 degrees per day; and 
the line cf apsides moves in the orbit plane about 0.85 degrees per 
day, therefore, the errors are small. 

3. The Vehicle Velocity Vector r^g 

The same assumptions are made for the velocity vector of the vehicle 

relative to the Earth as in the case of the vector r._ . The vector 

Vi!i 

• - 

r^ may be written as (Page 282 of Reference 8) 
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*Ve “ ^ “ p »f a ^ sin E + Q J 6 GM g U-e 2 ) cos E ) 

VE 


where. 


GMg * 5.986032 x 10 s k» 3 * sec” 2 - 

4. The Satellite-Sun Vector r^ 

The vector which points from the satellite to the Sun is 


r SV 3 r SE * r VE 

The vectors r_ and r„_, are given in sections 1 and 2 above. 


C. O RBITA L PARAMETERS AND SPACECRAFT ANGIES 

1 . True Anomaly v The true anomaly, 

is the angle between the vector from the center of the Earth to 
perigee and the vector from the center of the Earth to the satellite 


measured from the perifocus in the direction of motion. (See Figure 2). 

TYip following formulas were selected to calculate u (Page 3^1, 
Reference 9)s 


and 


cos v 


cos B - c 
1 - e cos £ 


sin v 


^L-e 2 sin E 
1-e cos E 
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2. Solar Array o 


The angle measured counter-clockwise from the vector to J 

vector when observed frcm the positive end of tne i^ vector is defined as 

the solar array angle, r . (See Figure l). Consistent with this notation 

P 

o o 

the planned range for the variation of a for OGO is 90 £ £ 270 . 

TP *P 

(Page 20, Beference 4). 6 is the angle between the hector and the k^ 

vector measured counter-clockwise from k^ when observed from the positive 

— o o 

end of the i b vector. It follows that when 9 goes from 90 to 270 , 6 

0 © F 

goes from 0 to 180 . Therefore, 9 = 5 + $0° and cos 6 = T • k. . 

p P o 


Therefore, tpp may be found from the formula 


-1 


c °s U * kj + 90 


where 


0° s cos" 1 (J * kjj) * 180° 


3. OPEP Angle $ e 

This angle is measured counter-clockwise from the vector i^ to the 

vector i e when observed frcm the positive end of the vector k^. (See 

Figure l). The OPEP control system is designed such that the range for 

o o 

the OPEP angle is ± 225 $ centered about = -90 . (Page kk, 

Beference 4). The analysis is simplified by assuming that the range 

_ o o 

for the OPEP angle is ± loO centered about t e 52 *90 « 





Define C 


It follows that cos If = i * i • 

e b e 

By geometrical considerations the quadrant of Y e 

the signs of cos V and C. 

e 


* * T e> ‘ S, : 

may he determined fran 


For cos i' e > 


0. Y is 

e . 


in the 1st or 4th quadrant. 


For cos Y 

e 

For C 


For C 


< 0, Y g is in the 2nd or 3rd quadrant. 


> 0, Y is in the 1st or 2nd quadrant. 


< 0, y is in the 3rd or 4th quadrant. 


4. Angle Between the Line of Apsides and the Ecliptic Plane Y 

This is the minimum angle between the vector pointing from the center of 

o 6 

the Earth to perigee and the ecliptic plane. This angle varies from 0 to 90 . 
Frcaa geometrical considerations it follows that: 

sin Y = P • k’ 

where 

k’ = -J sin « + k cos e 

+ IPj + 

? • k* = cos s - ~y E — s 
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5 


• Flight Path Angle y 

This is the min Irm a angle between the satellite's velocity vector ryg 

and the plane nor mal to r . The following formula was used to calculate 

VE 

V (Page 3h2, Reference 9)1 


tan y 


e sin y 
1 + e cos v 


6. OPEP - Velocity Vector Angle 6 

This is the angle between the vectors i and r._ . 

e Vi u 

e , cos- 1 

I *ve| 

where i^ is defined on page 12. 



7 . Angle Between the Projection of the Line of Apsides on the Ecliptic 
Plane and -the Earth-Sun Line 5, 

This is the angle between the projection of the -P vector (-P points from 
the center of the Earth to apogee) on the ecliptic plane and the vector 
which points froa the center of the Earth to the Sun. 


I = cos 


- 1 | -k' x(Pxk') ’ "SE 1 

I | k' x (P x k’ ) \ J 
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8. The Distance From the Center of the Earth to the Satellite, r^j. 

The following formula was used to calculate r^j, (Page 336, 
Reference 9); 

_ = a (1-e 2 ) 

VE 1 + e cos v 

D. HEAT INPUTS TO AN ARBITRARILY ORIENTED SURFACE OF EGO 
1. Earth Knitted Radiant Heat 

The treatment considers the Earth as a black body radiating 
uniformly at T = 250° K (Reference lo) . The emittance A is given 

by A = a T 4 , vhere a is the Stephan-BoTtzman constant and T is 
the absolute temperature for a black body. Therefore the average 
quantity of power radiating into space from the surface of the Earth 
is A * 0.02215 watts per square centimeter of the Earth’s surface. 

Some portion of the integrated heat flux radiated from the Earth will 
be incident upon a surface of the satellite whose unit normal is n 
(See Figure 3). The heat flux, q^, in integral form (See Figure 3 ) 
is: 

^ « A J fg (TU \ l , p) dA 
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where 


cos T| cos 


if cos T, > and cos u > 0. 


f E (T), p,, p) = 0 if cos ^ ^ 0 or cos p, « 0 


= An elemental area of the Earth’s surface. 

= The unit outward-directed normal of the surface in 
ijuestion. 

= The vector which points from the Earth’s elemental 
area dA to the satellite. 

= The angle between the outward-directed normal of dA 
and the vector p . If cos Tj s 0 then the satellite 
cannot see the area dA , hence f E s O. 

= The angle between the outward-directed normal of the 
surface, n and the vector -p which points from the 
surface to the area dA * If cos s 0, then the 
surface with normal n cannot see the area dA # hence 
f E *°. 


Define the Earth emitted radiation factor F-g by 


E — 


The values of Fj, have been evaluated in closed form 
(References 10 and ll). 

Depending upon the range of X (where X i s the angle between 
-r and n) there are three different formulas required for calculating 
Fg * They are as follows: 

a. If the surface of the satellite with outward normal n cannot see 
the Earth (see Figure **) 

. TT 

X * 2 + » m 

then F_ = 0 
E 

where 5 = sin** 1 ^ = the maximum look angle for which the 

2ft ll 

satellite can see the surface of the Earth* 


H = n , where r—, is the distance from the center of the Earth 

Yi* 

to the vehicle and R^ is the equatorial radius of the Earth. 


b. If 


~ $ < X < 


m 


a + 


a 
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V - l _ sin” 1 I 

r (a 2 -!)* 1 + 


£ ^ t. 1 

[_ fi sin X J 

' • 

1 i - 1 

^^2 C0S A cos 

r -(i^-i)^ cot x] 



£ l-H 2 cos 2 X 1^ 



For this ease, the surface is so oriented that its plane falls 
within the tangent cone (tne cone defining the limits of the visible 
area of the Earth) so that the surface with unit normal n receives 
radiation from oily a portion of the total area visible from the 
satellite. (See Figure 5)» 

e. If 

0 * X * | 

„ _ cos X 

F E ~ ~ 

For this case the surface is oriented such that it can 6ee 
entire spherical cap which is defined by the tangent cone. 

(See- Figure 6). 



2. Direct Solar Radiant Heat, — The heat flux incident upon the surface of 
the satellite whose unit normal, is n, which arrives directly from the Sun. Tne 


fraction of the solar constant S which is incident upon this surface of the 


F 


S 


n ' rsv 

‘PsvT 




satellite is 


S 



3* Reflected Solar 


Radiant Heat Qn -- Tne heat fl'ox incident upon the surface 
of the satellite vhose unit normal is n } which originates at the Sun and 
reflects from the Earth to the satellite. In integral form (see Figure 3). 

% = S c J f a P) dA 

cos T t cos u cos £ 

'• p 2 
if 

cos T\ > 0, cos y, > 0 and cos a > 0 

and 

\ (i), v, e, p) = 0 

if cos "H s 0 or cos y £ 0 or cos S $0 

c = the albedo = the fraction of the solar constant which is 
reflected away from the Earth, (o = .$k for this study. Reference 12). 

S = The angle between the vector u qT , and the outward-directed normal 
of the Earth's elemental surface area &A (See Figure 3 )• 

If cos 3^0, then that portion of the Earth's surface dA is 
not lighted by the Sun, f n = 0. 


vnere 

^ (^1* y* 3# p) = 


- 26 - 


A closed fora solution could net be found for this integral. Its 
value was approximated by assuming an average value for cos p. It was 
assumed that cos 5 (average) = cos 6g (See Figure 3) where 


®S = The angle between the unit vector u and the vector r, 

SE 


YE 


This reduces the integral to the same form as <jg. (see page 22). 
Hence, the approximate expression for a becomes 


\ = & S cos e s Pg 


where the closed form solution for F-, is presented on page 23. 


Define 

F c * 3 S F * C06 ®S 


except when 


cos 6g s 0 


then 

P * 0 
a 


because 

cannot be a negative quantity. 
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Figure 7 presents a comparison between q^ (approximate) ana q^ 
(actual) as a function of altitude with \ and 6 S as parameters. 

The values of q^ (approximate) were calculated using the simplified 
formula presented on page 27* The values of q^ (actual) were 
determined by numerical evaluation of the integral on page 26 of this 
report for 5 C = 180°. For this value of { c the difference between 
(approximate) and q^ (actual) is maximumized; the difference was 
calculated to be less than 2 m v/ cm*' as shown on Figure 7* The angle 
5 c is the angle between the plane formed by the vectors ryg and 
and the plane formed by the vectors end n (See Figure 8). The 
value of (actual) is a function of r^, X, 65 and $ c ; whereas q a 
(approximate) is a function of r^, X, 9g, but not a function of ? c . 
Figure 7 shows that, the maximum difference between the two methods 
occurs at the lower altitudes. The difference becomes insignificantly 
sm a ll at altitudes greater than 20,000 1¥.M. The formula for ^ 
(approximate) breaks down for 8g = 90°, but q^ (actual) becomes 
insignificantly small and can be approximated to zero. The q G 
approximation is adequate when compared to the uncertainties in the 
albedo, e.g. , the albedo varies between 0.3 and 0.6 depending on the 
^ latitude and season (Reference 13). The maximum difference between 
q c (approximate) and (actual) is only 1.3^ of the solar constant. 


COMPUTER PROGRAMS 

A. The Launch Window Program . 

The basic Launch Window Program was developed by P. Musen with 
the assistance of A. Smith (both of the Theoretical Division; Goddard 
Space Flight Center). R. Devaney, also of Goddard, made the changes for 
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its adaptation to the 0G0 problems. A description of the basic 
Launch Window Program, including a bibliography, is presented in 
Reference 1. 

The input quantities are the orbit elements a^, e Q , i Q , end 
^ at the initial tine T * The program calculates a, e, i, Q, ana 3 

O Q 

for time t, taking into account the influence of the gravitational 
forces due to the Moon, the Sun ana the Earth, including the 2nd, 3rd 
and bth zonal harmonics. The time per orbit that the satellite spends in 
the Earth’s shadow is also computed by the program. The program does 
not calculate the position of the spacecraft in its orbit. 

The orientation of the orbit in space is defined by the unit 
vectors P, Q and R. These are given in terms of i, 0, and a) in 
Reference 8. 

3. The Interplanetary Trajectory Program (IPP) 

This program (Reference lb) was used, to check the behavior of 
the orbital elements as a function of time. (Reference l). The XPP 
can be used to calculate the position and velocity of the satellite 
relative to the Earth as a function of time when given an initial set 
of injection conditions. The launch window program does not have this 
capability. The launch window program calculates the spacecraft 
angular data presented in this report, to a sufficiently high degree of 
accuracy that it may be used with confidence in most cases. All of the 
angle and heat input calculations presented in this report will be 
programmed into the IPP (with its accurate position calculation) and will 
be available for post flight analysis of the EGO (S-b9) satellite. 
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VI. IgSUUTS AMD DISCUSSION 

A. Orbital Parameters and Spacecraft Angles 

The results oi orbit parameter studies for EGO are shown graphically 
cm Figures 9 through l6. These studies are based on a launch date of 
November 6, 1963 at 16.0, 20.0 and 2k. 0 hours U.T. The elements of the 
initial orbits are: 


a = 

9.731167 

Earth radii 

e = 

.8929016 


i = 

30°. 807 


m - 

-45*°596 


p a = 

75-°27 

for 16.0 hours U.T. 

a> = 

135. 0 43 

for 20.0 hours U.T. 

°3 = 

195. °59 

for 24.0 hours U.T. 


These, elements are consistent with the initial conditions of 
Part XU on page 13. 

Due to the large semi-major axis (a = 9 *73317 Earth Radii), the 
period for EGO is 42.75 hours. The orbit has great variations in the 
speed of the spacecraft throughout its orbit, especially near perigee, 
as can be seen in the rapid variations in true anomaly., solar array 
angle, OEEP angle, flight path angle, OFEP-velocity vector angle, and 
in the distance from the center of the Earth to the satellite; shown 
on Figures 9, 10, 11, 13, 14 and l6, respectively. 
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1- Trae anomaly, p 

Figure 9 shews the true anomaly v ; ition with time for one 
orbital period. The confutation was done with two-body 
equations and the orbital elements at- the time of injection 
were used for this computation. This is a typical curve and 
its shape will change slightly throughout the year as the orbit 
changes due to perturbations. 

2. Solar array angle, 

Figure 10 shows the history of the solar array angle for # 

one orbital period. The angle history varies as a function 
of injection time into orbit due to the motion of the Earth and 
Sun. Figure 10 shows the angle histories for three different 
injection times. 

J 

3. OPEP angle, $ e . 

The OPEP angle variation with time is shown on Figure 11 and 
like the solar array angle its history is a function of 
injection time (a curve is shown for three injection times). 

h. Angle between the line of apsides and the ecl iptic plane, $ 

Figure 12 shows the variation of angle between the line of 
apsides and the ecliptic plane for one year for three different 
injection times of a given launch day. The absolute value of 
the expression for sin -J is chosen so that the angle shown on 

Figure 12 will always lie between 0° to 90°. flench, in this study 
the value of the angle $ is shown without regard to quadrant. 
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• Flight path angle, "Y 

The flight path angle varies periodically with each canplete 
orbit and changes only as the shape of the orbit changes. 

Figure 13 shows the flight path angle for one orbital period 
measured from perigee. The Launch Window Program has the 
capability of generating this angle at any time. The history of 
the angle for any orbit during the first year will not differ 
greatly from Fig. 13* 

6. OPEP - Ve locity Vector Angle, 

The variation of the angle between the Xg axis of orbit plane 
experiment package (OPEP) and the satellite* s velocity vector will 
be periodic with each orbital period. Figure 14 shows the angle 
as a function of time. This is a typical curve and will c han ge 
only slightly throughout the year. The shape of the curve is 
independent of launch time. Since the absolute value of the 
expression for 9 g is used, the angle is shown on the figure without 
regard to the quadrant (See Part XV, Section C-6). 

7 . Angle Between the Projection of the Line of Apsides on the 
Ecliptic Plane and the Earth-Sun Line, c, 

Since the motion of the perigee point (to) and the precession of the 
line of nodes (fl) always maintain the same sign the projection of 
the line of apsides on the ecliptic plane will continually revolve. 
Consequently, the angle between the projection of the line of 
apsides cm the ecliptic plane and the Earth-Sun line will continu a lly 
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vary from zero to 2-n radians. It facilitates computation to 
show this angle varying only through l 80 °, I. e. , the smaller of the 
- angles between the projection of the line of apsides on the 
ecliptic plane and the Earth-Sun line. (See Figure 15)* The 
curves an Figure 15 show the periodic angle variation for three 
different injection times of a given launch day. The different 
curves as a function of injection time are due to the change in 
position of the injection point with respect to the inertial frame 
because of the Earth's rotation. 

8. The distance from the center of the Earth to the satellite, r^j. 

Since the orbit for the EGO satellite is non-circular the distance 
between the spacecraft and the center of the Earth will vary 
periodically with each orbit. (See Figure 16). The ch ange in the 
periodicity and amplitude will vary only slightly during the 
first year of life of the satellite. 

B. Heat Inputs to EGO 

Figures 17 through 25 present the heat inputs incident upon the 
different faces (Figure l) of the satellite as a function of time from 
perigee for the first orbital period (42.75 hours). The data in these 
figures are based on a time of injection into orbit of November 6, 19^3 « 
at 24.0 hours U.T. These curves are intended to demonstrate the capability 
to calculate heat input s v rathe r than to provide design data. The Earth 


emitted and solar reflected beat inputs axe negligibly small when the 
satellite is "beyond an altitude of 20,000 nautical Giles or correspondingly 
2.63 hours flight tine from perigee. 

1. Main box 

a. Earth emitted 

Figure 17 presents the Earth emitted heat inputs (<3g) to 
the main box as a function of time perigee. The + z^ face has the 

largest Earth emitted heat input since it always faces the Earth. 
Conversely, the * z- Q face has a zero Earth emitted heat input. The 
+ x^, “ x b> + y b , end -y^ faces all have equal Earth emitted heat 
inputs due to geometrical symmetry with respect to the Earth. The 
curve shown is independent of launch time. 

b. Direct solar 

Figure 18 presents the direct solar heat inputs («lg) to the 
main box as a function of time from perigee. This is a typical curve. 

The solar heat inputs to the + x^, - and + y^ face are always zero 
(see page 12 of this report). 

Tue maximum solar input to any other face is 137 milliwatts/ 

p 

cm which would occur when the Sun’s rays are normal to a particular 
surface. 

c. Reflected solar 

Figure 19 presents the reflected solar heat inputs <v> to 
the main box as a function of time from perigee. The + z fa face has the 
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largest reflected solar heat input since it faces the Earth. The . 

- ^ face has zero reflected solar heat input. The + x^, - x^, + y b > 
_ faces all have equal reflected solar inputs since it is assumed 
that q a (appro ximat e) is independent of t c . For this l au n ch date 
(November 6, 1963, 24.0 hours U.T. ) the angle between r^ and r^ was 
= 63 . 3 ° at perigee and the solar reflected heats are q^ - 6.78 
millivatts/cm 2 da the + - x b , + y b and - y b faces; and 19-32 

milliwatt s/cm^ on the + z b face. 


2. Solar Array 

a. Earth emitted 

Figure 20 presents the Earth emitted heat inputs (<^) to the 
solar array as a function cf time from perigee. The + and - x^ faces 
of the solar array have the same Earth emitted heat inputs as the 
+ ^ and - face of the main boot (See Figure l). The Earth emitted 
heat inputs (<^,) to the + y^> “ + z p and " z p faces 6Xe 6150 

included on this figure. 


b. Direct solar 

Figure 21 present- :.e direct solar heat inputs (qg) to the 
solar array as a function of time from perigee. Since the + y^ fece 
of the solar array contains the solar cells and looks directly at the 
Sun (See page 11 of this report) the heat into this face is the full 
value of the solar constant S = 137 • 4 mil ll watts/ cm » The direct 
solar heat inputs to the other faces of the array - y^> + x^, - x p , + z^ 

* Z p are all zero* 
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c. Reflected solar 


Figure 22 presents the reflected solar heat Inputs (q ) to 

c 

the solar array as a function of time from perigee. The + x^ and - x^ 
faces of the solar array have the same reflected solar heat inputs 
as the + Xjj and - face of the main box. The reflected solar heat 
inputs to the + 7 p , - y p , + * p “* * *p •“ also *“ 

this figure. 


3 . Orbit Plane Experiment Package (OHSP) 
a. Earth emitted 

Figure 23 presents the Earth emitted heat inputs (q^,) to 

the OPEP as a function of time from perigee. The + z g face has the 

sane Earth emitted heat inputs as the + z^ face of the main box since 

both faces are oriented the same way. The - z face has a zero 

e 

Earth emitted heat input since it never faces the Earth. The + x g , 

- x^, + y g , - y g faces all have equal Earth emitted heat inputs due 
to geometrical symmetry with respect to the Earth. 


b. Direct solar 

Figure 2h presents the direct solar heat inputs ( q^) to the 
OPEP as a function of time from perigee. The solar heat inputs to the 
+ z and - z faces are the seme as to the + z^ and - z^ faces of the 
main box, respectively. The solar inputs to the + y^ and - y^ faces 
are nearly constant over one orbital period since these faces are 
parallel to the orbit plane and since the angle "between the Earth-Sun 
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line and the plane of the orbit is nearly constant over one orbital 
period. The direct solar heat inputs to the + x^ and -x e faces are 
also included on the figure. 

c. Reflected solar 

Figure 25 presents the reflected solar heat Inputs (q^) to 

the OEBP as a function of time from perigee. The + face has the 

same reflected solar heat inprut as the + face of the main box since 

both faces are oriented the same vay. The - z face has .zero solar 

e 

reflected heat input, since it never faces the Earth. The + x , - x> 

e e 

+ y , - y g all have equal solar reflected heat inputs due to geometric 
symmetry with respect to the Earth. 
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VII. CONCLUSIONS 


The methods of analyzing the orbital parameters, the spacecraft angles 
and the spacecraft heat inputs presented in this document were developed 
primarily to do a pre -launch study of the Eccentric Orbiting Geophysical 
Observatory, S-h9 (EGO), but the methods are general and may be applied to 
any near Earth satellite. These methods especially have application for studies 
on future Orbiting Geophysical Observatories (e.g., FOGO). 

The program developed as a result of this study is a modified Halphen 
method program. The basic program was conceived by P. A. Musen and 
subsequently modified for use as a launch window program for EGO. Flexibility 
in the choice of control laws makes the program appealing for general use. 

A secondary "but important advantage of the program is that a year study may 
be computed in about three minutes on -the IBM 709k. 

The closed form solution of the solar reflected heat input developed in 
this document, while being an approximation, was found to be. quite accurate. 

The maximum error resulting from the use of the closed form solution is only 
1. 3J& of the solar constant, and it is not restricted to the problem treated. 

The program has been named the "Launch Window Program". It has been 
designed for use on IBM 709 k computer but can be easily adapted for other 
computing systems. 


- 38 - 


VIII. references 


1. Montgomery, H. E. , EGO Launch Window Study, NASA Report 
X643-62-225, Goddard Space Flight Center, Greenhelt, Maryland. 

2. Hemenover, A. P. , Nominal Trajectory and Impact Data of EOGO 
Mission, Lockheed Missiles and Space Company Report No. A-306775/ 
91-21, dated November 9> 19^2 , Sunnyvale, California. 

3* Specification for the Orbiting Geophysical Observatory, S-49/50 
Revision 1, dated August 19 62, Goddard Space Flight Center, 
Greenbelt, Maryland. 

4. Otten, D. D. , OGO Attitude Control Subsystem Description, 

Logic, and Specifications, Report Number 2313-0004-HV-000, 

4 December 1961, Space Technology Laboratories, Inc., Redondo 
Beach, California. 

5. Moulton, F. R. , Celestial Ifechanics, Macmillan, New York, 1914. 

6. Clarke, Victor C., Constants and Related Data Used in Trajectory 
Calculations at the Jet Propulsion Laboratory, Technical Report 
Number 32-273. May 1, 1962, Jet Propulsion Laboratory, California 
Institute of Technology, Pasadena, California. 

7« The American Ephemeris and Nautical Almanac, i960. United 
States Government Printing Office, Washington, D. C. 

8. Baker, Robert M. L. and Makemson, Maud W. , An Introduction to 
Astrodynaraics , Academic Press, New York. 

9. Kraft, Ehricke A., Space Flight, Volume 1, D. Van Nostrand . 
Company, Inc. 

+ 

Cunningham, F. G. , Power Input to a Small Flat Plate From a 
Diffusely Radiating Sphere, With Application to Earth Satellites, 
NASA TN D-710. August, 1961, Goddard Space Fligit Center, 
Greenbelt, Maryland. 


• 39 - 


10 . 


11. Stevenson, J. A. and Grafton, J. C., Radiation Heat Transfer 
Analysis for Space Vehicles, ASP Technical Report 6l-ng, 
Space and Information Systems Division, North American 
Aviation, Inc. 


12. Allen, C. W., Astrophysical Quantities, University of London, 
The Athlone Press, 1955/ Published in the U.S.A. by Essential 
Books, Fairlavn, Nev Jersey. 


13. Johnson, John C., Physical Meteorology, John Wiley and Sons, 
Inc., Nev York, 1954. 


14. Pines, S. , and Wolf, H., Interplanetary Trajectory by the 
Encke Method Programmed for the IEM 7090, Report No. RAC - 
656-451» 12 December, i960. Research and Development Division 
Republic Aviation Corporation, Farming dale. Long Island, 

Hev York. 


- 40 - 





CIRCLE CENTERED AT 
CENTER OF ELLIPSE 




CASE 1. SATELLITE SURFACE WITH UNIT NORMALn” 
CANNOT SEE THE EARTH. 




CASE 2. SATELLITE WITH UNIT NORMAL n CAN SEE A PORTION 
OF THE EARTH’S SURFACE WITHIN THE TANGENT CONE. 



Fig. 5 Geometry for Earth 


d Heat Inputs 


CASE 3. SATELLITE SURFACE WITH UNIT NORMAL ttCAtf SEE 
THE ENTIRE SURFACE WITHIN THE TANGENT CONE. 



Fig. 6 Geometry for Earth Emitted H 


uts 


REFLECTED SOLAR HEAT INPUT, q<* mw/em* 





Fig. 8 Geometry (or Reflected Solar Heat Inputs 
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Figure 9. True anomaly 



DATE OF INJECTION INTO ORBIT NOV. 6, 1963 
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Figure 10. Solar array angle 
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Figure 12. Apsides-ecliptic angle 





TIME OF INJECTION INTO ORBIT: 
NOV. 6, 1963, 16.0 HOURS U.T. 
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Figure 13. Flight path angle. 




TIME OF INJECTION INTO ORBIT: 
NOV. 6, 1963, 16.0 HOURS U.T. 
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Figure 14. Opep-velocity vector angle. 
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Figure 15. Angle between the projection of the line of 

apsides on the ecliptic plane and the 
earth-sun line. 
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